Experimental determination of the evolution of the Bjorken integral at low Q 2 
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We extract the Bjorken integral Y\~ n in the range 0.17 < Q 2 < 1.10 GeV 2 from inclusive scattering 
of polarized electrons by polarized protons, deuterons and 3 He, for the region in which the integral 
is dominated by nucleon resonances. These data bridge the domains of the hadronic and partonic 
descriptions of the nucleon. In combination with earlier measurements at higher Q 2 , we extract the 
non-singlet twist-4 matrix element /2. 
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For almost 50 years experimental and theoretical re- 
search efforts in hadronic physics have sought to under- 
stand the structure of the nucleon. With the development 
of Quantum Chromodynamics (QCD), these studies have 
focused on obtaining an accurate description of nucleon 
structure in terms of fundamental quark and gluon de- 
grees of freedom. A powerful tool has been deep inelastic 
lepton scattering from nucleons and nuclei, and the as- 
sociated theoretical machinery of the operator product 
expansion (OPE), which allows the interpretation of the 
measured structure functions in terms of parton momen- 
tum and spin distribution functions. 

Experiments using polarized beams and targets have 
played a critical role in testing the application of QCD 
to nucleon structure 0. The Bjorken sum rule 0], which 
relates the first moment of polarized deep inelastic struc- 
ture functions to nucleon ground state properties, has 
been an important part of these studies. At infinite four- 
momentum transfer squared, Q 2 , the sum rule reads 
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where g\ and g™ are the spin-dependent proton and neu- 
tron structure functions, respectively. Here, qa is the 
nucleon axial charge and x = Q 2 /2Mv, with v the en- 
ergy transfer and M the nucleon mass. The sum rule has 
been verified experimentally to better than 10% 0,0,0 
Because the Bjorken sum rule relates differences of 
proton and neutron structure function moments, and 
r™, only flavor non-singlet quark operators appear in the 
OPE. Another simplification arises at low Q 2 where the 



resonance contributions to the proton and neutron, in 
particular that of the A(1232) resonance, partly cancel. 
This cancellation simplifies calculations based on chiral 
perturbation theory (%PT), and may extend the Q 2 range 
of their applicability. The gap between the domains of 
validity for perturbative QCD (pQCD) and %PT might 
even be bridged, enabling for the first time a fundamen- 
tal theoretical description of nucleon structure from large 
to small scales 0. The Bjorken sum rule is therefore 
relevant for understanding the transition from pQCD to 
nonperturbative QCD. 

In this Letter we report on a determination of the 
Bjorken integral using data obtained at the Thomas Jef- 
ferson National Accelerator Facility (Jefferson Lab) over 
the Q 2 range of 0.17-1.10 GeV 2 . Combined with higher 
Q 2 data from earlier experiments, we analyze the data 
using the OPE and extract the l/Q 2 higher twist correc- 
tions to the integral at intermediate values of Q 2 . 

The data were obtained in three different experiments 
using polarized electrons on polarized proton 0, deu- 
terium and 3 He targets. To analyze the 
scattered electrons, the proton and deuteron experi- 
ments used the CEBAF Large Acceptance Spectrometer 
(CLAS) in Hall B 0, while the 3 He experiment used 
the two High Resolution Spectrometers in Hall A (l2^ . 

The individual measurements of the proton, neutron 
and deuteron integrals Y\ ,n ' d have been reported else- 
where E3 ■ To form the isovector combination 
T\~ n we subtract from the experimental values of g\ the 
values of extracted from the 3 He or the deuteron mea- 
surements. However, in order to combine these data, the 
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TABLE I: Inelastic contributions to the Bjorken sum. The 
second and third columns give the sum and its uncertainty for 
W < 2 GeV. The fourth to sixth columns give the total sum 
and its uncertainties. The last column indicates the origin of 
the neutron information (from 3 He or from deuteron data). 
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Q 2 values at which and g\ were obtained must coin- 
cide. We chose to re-analyze the 3 He data at six values 
of Q 2 which match the ones of the prot on data and differ 
from the values reported in Refs. [fSLllCj. For the deuteron 
measurement, given the larger uncertainties, we simply 
interpolated the proton data points to match the four Q 2 
points of the deuteron data. The additional systematic 
uncertainty from the interpolation is negligible. 

The three experiments 0, H El E3 have measured g\ 
up to an invariant mass W — 2 GeV. The unmeasured 
contributions to the proton and neutron integrals, corre- 
sponding to the low- a; domain, need to be consistently ac- 
counted for. In the current analysis, the fit from Bianchi 
and Thomas |l3j was used to estimate the low- a; contri- 
bution to the moments up to W 2 = 1000 GeV 2 . The 
uncertainty on this contribution was evaluated by taking 
the quadratic sum of the differences induced by indepen- 
dently varying each parameter of the fit within the range 
of values given in 13]. The contribution beyond 1000 
GeV 2 was determined using a Regge parametrization 
constrained so that the Bjorken integral at Q 2 — 5 GeV 2 , 
from the world data complemented by the Bianchi and 
Thomas fit and our Regge parametrization, agrees with 
the sum rule. The systematic uncertainties from the neu- 
tron and proton data have been added in quadrature. 
The moment T™ was extracted from 3 He or deuterium 
data usin g th e formalism of nucleon effective polariza- 
tions [313. The resulting T\~ n is shown in Fig. □ 
by the filled symbols, with the values given in Table I. 
Note that only the inelastic contributions are included 
in rf" n . Data from the SLAC E143 experiment [lfj are 
also shown (open circles) for comparison. 

The data are compared with theoretical calculations 
based on xPT and with phenomenological models. At 
Q 2 = 0, the slope of the Bjorken integral is constrained 
by the Gerasimov-Drell-Hearn (GDH) sum rule 0, • 
The Soffer-Teryaev model agrees only with the low- 
Q 2 data. The overestimate at larger Q 2 was traced back 



FIG. 1: Inelastic contribution to the Bjorken sum. The trian- 
gles (squares) represent the results with the neutron extracted 
from 3 He (deuteron) data, with the horizontal bands cover- 
ing the 0.17-0.99 GeV 2 and 0.34-1.1 GeV 2 ranges the corre- 
sponding systematic uncertainties. The E143 data are shown 
for comparison. The gray band represents the leading-twist 
(LT) contribution calculated to third order in a 3 . The curves 
correspond to theoretical calculations (see text). 



to the QCD radiative corrections and has now been cor- 
rected [2t|. The Burkert-Ioffe model [2lJ agrees well 
with the data over the full range covered in Fig. 1. This 
may indicate that vector meson dominance gives a rea- 
sonable picture of the parton-hadron transition. At low 
Q 2 attempts have also been made to calculate T^ - " us- 
ing xPT [i^, Hll- The calculation done in the heavy 
baryon approximation |23| seems to agree better with 
the data. At higher Q 2 the data are compared with the 
leading- twist calculation (gray band in Fig. 1), which cor- 
responds to incoherent scattering from individual quarks. 
In pQCD, gluon radiation causes scaling violations in 
structure functions, and introduces an a s dependence 
on the right hand side of Eq. JIJ. At leading- twist, the 
pQCD result at third order in a s (in the MS scheme) is 
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The gray band in Fig. 1 represents the uncertainty in 
™ due to the uncertainty in a s . There is reasonable 
agreement between the leading-twist prediction and the 
data. Their difference is related to higher twist effects 
that should become important at low Q 2 . In particular, 
application of the OPE to moments of structure func- 
tions requires the expansion of the total moment rather 
than the inelastic moment as in Fig. 1. While the elastic 
contribution is negligible at high Q 2 , it dominates at low 
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Q 2 . Fig. 2 shows the total moment, including the elastic 
contribution, calculated from the form factor parameter- 
izations from Ref. |24|. In addition to the JLab data, 
we also plot data at higher Q 2 from the SLAC E143 |25j 
and E155 0, DESY HERMES HJ and CERN SMC H 
experiments. For consistency, the low-s contributions, 
outside of the measured regions, have been re-evaluated 
using the same procedure as described earlier. 

The OPE analysis allows one to expand the total mo- 
ment r^ - ' 1 in powers of 1/Q 2 : 

oo p— n 

where the leading twist i = 1 coefficient is given in 
Eq. ©. The coefficients &7 n for % > 1 represent matrix 
elements of higher twist operators. The matrix elements 
contain information on the long range, nonperturbative 
interactions or correlations between partons. In particu- 
lar, the 1/Q 2 correction term is [27L Esj 

M 2 

tfT = ^- {at n + 44"" + 4/F") , (4) 

where a^~ n is the target mass correction given by the x 2 - 
weighted moment of the leading-twist g± structure func- 
tion, and g?2 _ ™ is a twist-3 matrix element given by 

<P~ n = [ dx x 2 (2«?r n + 3^~") ■ (5) 
"Jo 

The twist-4 contribution, f%~ , given by a mixed quark- 
gluon operator, is related to the color electric and mag- 
netic polarizabilities of the nucleon 

0,113. 

The a,2~ n correction, which is twist-2, is calculated us- 
ing the fit of polarized parton distributions from Ref. [3lJ • 
The d,2~ n matrix element is obtained from J^]. With 
these inputs, the data on T\~ n in Fig. 2 can be used 
to extract /f - ™. As an additional parameter in the fit, 
we include the 1/Q 4 coefficient n%~ n ■ For the leading 
twist contribution, we constrain the low- a; extrapolation 
by assuming the validity of the Bjorken sum rule for 
Q 2 > 5 GeV 2 . In fact, our low- a; extrapolation gives 
= 1.270 ± 0.045, which is very close to the empirical 
value gA = 1.267±0.004. The higher twist terms are then 
determined from the Q 2 < 5 GeV 2 data using our fitted 
value of gA ■ The point to point correlated uncertainty for 
the JLab data extracted from 3 He and hydrogen has been 
disentangled from the uncorrelated uncertainty, and only 
the latter is used in the fit. The correlated systematics 
are then propagated to the fit result, as is the uncertainty 
arising from a s . The data from the other experiments are 
treated as uncorrelated from point to point. 

It is not clear a priori over which Q 2 range the 1 /Q 2 ex- 
pansion should be valid. For instance, at Q 2 s» 0.7 GeV 2 
the elastic and leading twist contributions are of com- 
parable magnitude. Fitting over the range 0.8 < Q 2 < 
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FIG. 2: Total moment Fj"™, including the inelastic contribu- 
tion from Fig. 1 together with the elastic. The data extracted 
from 3 He (deuterium) together with proton data are indicated 
by the triangles (squares). The leading-twist (LT) contribu- 
tion is given by the gray band. The point to point correlated 
uncertainty for the data extracted from 3 He and proton is 
shown by the horizontal band. The error bars on the symbols 
represent the uncorrelated uncertainty. The data from other 
experiments are assumed to be uncorrelated. The fits to the 
total moment are indicated by the solid curves. 

10 GeV 2 gives jf" = -0.13 ± OJ5(micor)j:8;§f (cor), 
normalized at Q 2 = 1 GeV 2 , where the first and sec- 
ond errors are uncorrelated and correlated, respectively, 
and l4~ n /M i = 0.09 ± 0.06(uncor) ± 0.01(cor). The 
contribution to the total uncertainty from the elastic 
form factors is negligible. Starting the fit at a lower 
Q 2 , Q 2 — 0.66 GeV 2 , yields the more negative value 
jp-n _ _o.l8±0.05(uncor)to^5( cor )' anc ^ a larger value 
for i4~ n , fj%- n /M 4 = 0.12 ± 0.02(uncor) ± 0.01 (cor), 
with somewhat smaller errors. The results of the two fits 
are shown in Fig. 2, but are almost indistinguishable. 
At Q 2 = 1 GeV 2 , the 1/Q 4 contribution is M6~"/Q 4 - 
0.09 ± 0.02, which is of similar magnitude and of oppo- 
site sign to the 1/Q 2 term, ^ P f n /Q 2 ~ -0.06 ± 0.02, 
obtained by adding the extracted /f - ™ value to d^"™ 
and aP 2 ~ n in Eq. (4). This may explain why the leading 
twist description agrees well with the data down to sur- 
prisingly small values of Q 2 (~ 1 GeV 2 ), and could be a 
hint that quark-hadron duality might work well for p — n 
non-singlet quantities. 

These results also suggest that at these lower Q 2 val- 
ues the twist expansion may not converge very quickly, 
and that higher twist terms may be needed. Including a 
jUg n /Q 6 term, however, gives significantly larger uncer- 
tainties on the higher twist contributions, making them 
compatible with zero. Starting the fit at Q 2 = 0.47 GeV 2 , 
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for instance, gives 
+0.04(cor), 



-0.14 ± O.lO(uncor) ± 



/M 4 = 0.09±0.08(uncor)±°;^(cor) and 



- n /M e = 0.01 ± 0.03(uncor) ± 0.02(cor). 



Ms 

The results for f%~ can be compared to non- 
perturbative model predictions: f%~ n = —0.024 ± 0.012 
Q and H~ n = -0.032 ±0.051 (QCD sum rules), 
fP_ n = 0.028 H3 (MIT bag model) and /| _n = -0.081 
[25j (instanton model). The results can also be com- 
pared with values obtained in analyses of the proton j3|| 
and neutron data separately. Naively subtracting 
f£ from f% gives 0.01 ± 0.08, which is consistent within 
uncertainties with the above values for However, 
different Q 2 ranges were used in the proton and neutron 
fits, and different low- a; extrapolations implemented. 

The larger uncertainty on f% n from the Q 2 > 
0.8 GeV 2 analysis reflects the larger values of Q^in used 
here compared with that used in the neutron analy- 
sis [3]) . Fitting the neutron data from Q 2 = 1 GeV 2 



rather from Q — 0.5 GeV as in Ref. 



would in- 



crease the uncertainty on f£ appreciably, which, when 
combined with the proton data fitted over the same 
range, would be more compatible with the uncertainty 
from the present combined analysis. This issue could 
be ameliorated with better quality data at higher Q 2 
(Q 2 > 1 GeV 2 ). Data in this region on the proton 
and deuteron collected in Hall B at Jefferson Lab are 
presently being analysed. Plans for high-precision mea- 
surements of the proton and neutron structure functions 
at higher Q 2 are also included in the 12 GeV energy up- 
grade of Jefferson Lab |38j . 

To summarize, we have presented an extraction of the 
Bjorken sum in the 0.17 < Q 2 < 1.10 GeV 2 range. Be- 
ing a nonsinglet quantity, the Bjorken sum simplifies the 
theoretical analyses at both high Q 2 (using the OPE) 
and at low Q 2 (using xPT). It thus provides us with 
a unique opportunity to understand better the transi- 
tion from perturbative QCD to the confinement region. 
Combining with data at higher Q 2 , we have extracted 
the higher twist contributions to the sum. We find f% ~ n 
small and the total higher twist contribution, for twists 
lower than eight, compatible with zero. 
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